Observational Constraints on the Response of High‐Latitude Northern Forests to Warming

Since the 1960s, carbon cycling in the high‐latitude northern forest (HLNF) has experienced dramatic changes: Most of the forest is greening and net carbon uptake from the atmosphere has increased. During the same time period, the CO₂ seasonal cycle amplitude (SCA) has increased by ~50% or more. Disentangling complex processes that drive these changes has been challenging. In this study, we substitute spatial sensitivity to temperature for time to quantify the impact of temperature increase on gross primary production (GPP), total ecosystem respiration (TER), the fraction of Photosynthetic Active Radiation (fPAR), and the resulted contribution of these changes in amplifying the CO₂ SCA over the HLNF since 1960s. We use the spatial heterogeneity of GPP inferred from solar‐induced chlorophyll Fluorescence in combination with net ecosystem exchange (NEE) inferred from column CO₂ observations made between 2015 and 2017 from NASA's Orbiting Carbon Observatory‐2. We find that three quarters of the spatial variations in GPP can be explained by the spatial variation in the growing season mean temperature (GSMT). The long term hindcast captures both the magnitude and spatial variability of the trends in observed fPAR. We estimate that between 1960 and 2010, the increase in GSMT enhanced both GPP and the SCA of NEE by ~20%. The calculated enhancement of NEE due to increase in GSMT contributes 56–72% of the trend in the CO₂ SCA at high latitudes, much larger than simulations by most biogeochemical models

Global Analysis of Bioclimatic Controls on Ecosystem Productivity Using Satellite Observations of Solar-Induced Chlorophyll Fluorescence

Ecosystem productivity models rely on regional climatic information to estimate temperature and moisture constraints influencing plant growth. However, the productivity response to these environmental factors is uncertain at the global scale and has largely been defined using limited observations from sparse monitoring sites, including carbon flux towers. Recent studies have shown that satellite observations of Solar-Induced chlorophyll Fluorescence (SIF) are highly correlated with ecosystem Gross Primary Productivity (GPP). Here, we use a relatively long-term global SIF observational record from the Global Ozone Monitoring Experiment-2 (GOME-2) sensors to investigate the relationships between SIF, used as a proxy for GPP, and selected bio-climatic factors constraining plant growth at the global scale. We compared the satellite SIF retrievals with collocated GPP observations from a global network of tower carbon flux monitoring sites and surface meteorological data from model reanalysis, including soil moisture, Vapor Pressure Deficit (VPD), and minimum daily air temperature (Tmin). We found strong correspondence (R2 \u3e 80%) between SIF and GPP monthly climatologies for tower sites characterized by mixed, deciduous broadleaf, evergreen needleleaf forests, and croplands. For other land cover types including savanna, shrubland, and evergreen broadleaf forest, SIF showed significant but higher variability in correlations between sites. In order to analyze temperature and moisture related effects on ecosystem productivity, we divided SIF by photosynthetically active radiation (SIFp) and examined partial correlations between SIFp and the climatic factors across a global range of flux tower sites, and over broader regional and global extents. We found that productivity in arid ecosystems is more strongly controlled by soil water content to an extent that soil moisture explains a higher proportion of the seasonal cycle in productivity than VPD. At the global scale, ecosystem productivity is affected by joint climatic constraint factors so that VPD, Tmin, and soil moisture were significant (p \u3c 0.05) controls over 60, 59, and 35 percent of the global domain, respectively. Our study identifies and confirms dominant climate control factors influencing productivity at the global scale indicated from satellite SIF observations. The results are generally consistent with climate response characteristics indicated from sparse global tower observations, while providing more extensive coverage for verifying and refining global carbon and climate model assumptions and predictions

Towards a harmonized long‐term spaceborne record of far‐red solar induced fluorescence

Far‐red solar‐induced chlorophyll fluorescence (SIF) has been retrieved from multiple satellites with nearly continuous global coverage since 1996. Multiple official and research‐grade retrievals provide a means for cross validation across sensors and algorithms, but produces substantial variation across products due to differences in instrument characteristics and retrieval algorithm. The lack of a consistent, calibrated SIF data set hampers scientific interpretation of planetary photosynthesis. NASA's Orbiting Carbon Observatory 2 (OCO‐2) offers small sampling footprints, high data acquisition, and repeating spatially resolved targets at bioclimatically diverse locations, providing a unique benchmark for spaceborne sensors traceable to ground data. We leverage overlap between the longer running Global Ozone Monitoring Instrument version 2 (GOME‐2) SIF time series, and more recent state‐of‐the‐art OCO‐2 and TROPOspheric Monitoring Instrument (TROPOMI) data, in a first attempt to reconcile inconsistencies in the long‐term record. After screening and correcting for key instrument differences (time of day, wavelength, Sun‐sensor geometry, cloud effects, footprint area), we find that Global Ozone Monitoring Instrument version 2 and TROPOspheric Monitoring Instrument perform exceedingly well in capturing spatial, seasonal, and interannual variability across OCO‐2 targets. However, Global Ozone Monitoring Instrument version 2 retrieval methods differ by up to a factor of 2 in signal‐to‐noise and magnitude. Magnitude differences are largely attributed to retrieval window choice, with wider windows producing higher magnitudes. The assumed SIF spectral shape has negligible effect. Substantial research is needed to understand remaining sensitivities to atmospheric absorption and reflectance. We conclude that OCO‐2 and TROPOspheric Monitoring Instrument have opened up the possibility to produce a multidecadal SIF record with well‐characterized uncertainty and error quantification for overlapping instruments, enabling back‐calibration of previous instruments and production of a consistent, research‐grade, harmonized time series

Towards a harmonized long‐term spaceborne record of far‐red solar induced fluorescence

Far‐red solar‐induced chlorophyll fluorescence (SIF) has been retrieved from multiple satellites with nearly continuous global coverage since 1996. Multiple official and research‐grade retrievals provide a means for cross validation across sensors and algorithms, but produces substantial variation across products due to differences in instrument characteristics and retrieval algorithm. The lack of a consistent, calibrated SIF data set hampers scientific interpretation of planetary photosynthesis. NASA's Orbiting Carbon Observatory 2 (OCO‐2) offers small sampling footprints, high data acquisition, and repeating spatially resolved targets at bioclimatically diverse locations, providing a unique benchmark for spaceborne sensors traceable to ground data. We leverage overlap between the longer running Global Ozone Monitoring Instrument version 2 (GOME‐2) SIF time series, and more recent state‐of‐the‐art OCO‐2 and TROPOspheric Monitoring Instrument (TROPOMI) data, in a first attempt to reconcile inconsistencies in the long‐term record. After screening and correcting for key instrument differences (time of day, wavelength, Sun‐sensor geometry, cloud effects, footprint area), we find that Global Ozone Monitoring Instrument version 2 and TROPOspheric Monitoring Instrument perform exceedingly well in capturing spatial, seasonal, and interannual variability across OCO‐2 targets. However, Global Ozone Monitoring Instrument version 2 retrieval methods differ by up to a factor of 2 in signal‐to‐noise and magnitude. Magnitude differences are largely attributed to retrieval window choice, with wider windows producing higher magnitudes. The assumed SIF spectral shape has negligible effect. Substantial research is needed to understand remaining sensitivities to atmospheric absorption and reflectance. We conclude that OCO‐2 and TROPOspheric Monitoring Instrument have opened up the possibility to produce a multidecadal SIF record with well‐characterized uncertainty and error quantification for overlapping instruments, enabling back‐calibration of previous instruments and production of a consistent, research‐grade, harmonized time series

Optimal model complexity for terrestrial carbon cycle prediction

The terrestrial carbon cycle plays a critical role in modulating the interactions of climate with the Earth system, but different models often make vastly different predictions of its behavior. Efforts to reduce model uncertainty have commonly focused on model structure, namely by introducing additional processes and increasing structural complexity. However, the extent to which increased structural complexity can directly improve predictive skill is unclear. While adding processes may improve realism, the resulting models are often encumbered by a greater number of poorly determined or over-generalized parameters. To guide efficient model development, here we map the theoretical relationship between model complexity and predictive skill. To do so, we developed 16 structurally distinct carbon cycle models spanning an axis of complexity and incorporated them into a model–data fusion system. We calibrated each model at six globally distributed eddy covariance sites with long observation time series and under 42 data scenarios that resulted in different degrees of parameter uncertainty. For each combination of site, data scenario, and model, we then predicted net ecosystem exchange (NEE) and leaf area index (LAI) for validation against independent local site data. Though the maximum model complexity we evaluated is lower than most traditional terrestrial biosphere models, the complexity range we explored provides universal insight into the inter-relationship between structural uncertainty, parametric uncertainty, and model forecast skill. Specifically, increased complexity only improves forecast skill if parameters are adequately informed (e.g., when NEE observations are used for calibration). Otherwise, increased complexity can degrade skill and an intermediate-complexity model is optimal. This finding remains consistent regardless of whether NEE or LAI is predicted. Our COMPLexity EXperiment (COMPLEX) highlights the importance of robust observation-based parameterization for land surface modeling and suggests that data characterizing net carbon fluxes will be key to improving decadal predictions of high-dimensional terrestrial biosphere models.</p

Influence of ENSO and the NAO on terrestrial carbon uptake in the Texas-northern Mexico region

Climate extremes such as drought and heat waves can cause substantial reductions in terrestrial carbon uptake. Advancing projections of the carbon uptake response to future climate extremes depends on (1) identifying mechanistic links between the carbon cycle and atmospheric drivers, (2) detecting and attributing uptake changes, and (3) evaluating models of land response and atmospheric forcing. Here, we combine model simulations, remote sensing products, and ground observations to investigate the impact of climate variability on carbon uptake in the Texas‐northern Mexico region. Specifically, we (1) examine the relationship between drought, carbon uptake, and variability of El Niño–Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO) using the Joint UK Land‐Environment Simulator (JULES) biosphere simulations from 1950–2012, (2) quantify changes in carbon uptake during record drought conditions in 2011, and (3) evaluate JULES carbon uptake and soil moisture in 2011 using observations from remote sensing and a network of flux towers in the region. Long‐term simulations reveal systematic decreases in regional‐scale carbon uptake during negative phases of ENSO and NAO, including amplified reductions of gross primary production (GPP) (−0.42 ± 0.18 Pg C yr^(−1)) and net ecosystem production (NEP) (−0.14 ± 0.11 Pg C yr^(−1)) during strong La Niña years. The 2011 megadrought caused some of the largest declines of GPP (−0.50 Pg C yr^(−1)) and NEP (−0.23 Pg C yr^(−1)) in our simulations. In 2011, consistent declines were found in observations, including high correlation of GPP and surface soil moisture (r = 0.82 ± 0.23, p = 0.012) in remote sensing‐based products. These results suggest a large‐scale response of carbon uptake to ENSO and NAO, and highlight a need to improve model predictions of ENSO and NAO in order to improve predictions of future impacts on the carbon cycle and the associated feedbacks to climate change

Carbon Monitoring System Flux Net Biosphere Exchange 2020 (CMS-Flux NBE 2020)

Here we present a global and regionally resolved terrestrial net biosphere exchange (NBE) dataset with corresponding uncertainties between 2010-2018: Carbon Monitoring System Flux Net Biosphere Exchange 2020 (CMS-Flux NBE 2020). It is estimated using the NASA Carbon Monitoring System Flux (CMS-Flux) top-down flux inversion system that assimilates column CO2 observations from the Greenhouse Gases Observing Satellite (GOSAT) and NASA\u27s Observing Carbon Observatory 2 (OCO-2). The regional monthly fluxes are readily accessible as tabular files, and the gridded fluxes are available in NetCDF format. The fluxes and their uncertainties are evaluated by extensively comparing the posterior CO2 mole fractions with CO2 observations from aircraft and the NOAA marine boundary layer reference sites. We describe the characteristics of the dataset as the global total, regional climatological mean, and regional annual fluxes and seasonal cycles. We find that the global total fluxes of the dataset agree with atmospheric CO2 growth observed by the surface-observation network within uncertainty. Averaged between 2010 and 2018, the tropical regions range from close to neutral in tropical South America to a net source in Africa; these contrast with the extra-tropics, which are a net sink of 2.5±0.3 Gt C/year. The regional satellite-constrained NBE estimates provide a unique perspective for understanding the terrestrial biosphere carbon dynamics and monitoring changes in regional contributions to the changes of atmospheric CO2 growth rate. The gridded and regional aggregated dataset can be accessed at http://doi.org/10.25966/4v02-c391 (Liu et al., 2020).

Improved Constraints on Northern Extratropical CO2 Fluxes Obtained by Combining Surface-Based and Space-Based Atmospheric CO2 Measurements

© 2020. The Authors. Top-down estimates of CO2 fluxes are typically constrained by either surface-based or space-based CO2 observations. Both of these measurement types have spatial and temporal gaps in observational coverage that can lead to differences in inferred fluxes. Assimilating both surface-based and space-based measurements concurrently in a flux inversion framework improves observational coverage and reduces sampling related artifacts. This study examines the consistency of flux constraints provided by these different observations and the potential to combine them by performing a series of 6-year (2010–2015) CO2 flux inversions. Flux inversions are performed assimilating surface-based measurements from the in situ and flask network, measurements from the Total Carbon Column Observing Network (TCCON), and space-based measurements from the Greenhouse Gases Observing Satellite (GOSAT), or all three data sets combined. Combining the data sets results in more precise flux estimates for subcontinental regions relative to any of the data sets alone. Combining the data sets also improves the accuracy of the posterior fluxes, based on reduced root-mean-square differences between posterior flux-simulated CO2 and aircraft-based CO2 over midlatitude regions (0.33–0.56 ppm) in comparison to GOSAT (0.37–0.61 ppm), TCCON (0.50–0.68 ppm), or in situ and flask measurements (0.46–0.56 ppm) alone. These results suggest that surface-based and GOSAT measurements give complementary constraints on CO2 fluxes in the northern extratropics and can be combined in flux inversions to improve constraints on regional fluxes. This stands in contrast with many earlier attempts to combine these data sets and suggests that improvements in the NASA Atmospheric CO2 Observations from Space (ACOS) retrieval algorithm have significantly improved the consistency of space-based and surface-based flux constraints